Diffusion in the special theory of relativity 
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The Markovian diffusion theory is generalized within the framework of the special theory of 
relativity using a modification of the mathematical calculus of diffusion on Riemannian manifolds 
(with definite metric) to describe diffusion on Lorentzian manifolds with an indefinite metric. A 
generalized Langevin equation in the fiber space of position, velocity and orthonormal velocity frames 
is defined from which the generalized relativistic Kramers equation in the phase space in external 
force fields is derived. The obtained diffusion equation is invariant under Lorentz transformations 
and its stationary solution is given by the Jiittner distribution. Besides a non-stationary analytical 
solution is derived for the example of force-free relativistic diffusion. 

PACS numbers: 05.10.Gg, 03.30.+p, 04.20.-q 



1. INTRODUCTION 

The formulation of a consistent theory of Markovian 
diffusion within the framework of the relativity theory 
is a long-standing problem in physics. Over the years 
umpteen studies have been devoted to this issue (see 
e.g. Q - [13) with different and apparently irreconcil- 
able points of view. At least up to the knowledge of 
the author, a general accepted consistent solution of this 
problem is still outstanding. However, besides its fun- 
damental theoretical interest such theory is of particular 
importance in several applications such as in high energy 
collision experiments (see e.g. [13)) astrophysics (see e.g. 
[Hj]) and others. 

An alternative physical approach for the description 
of a relativistic gas in a heat bath is given by statistical 
thermodynamics and the Boltzmann equation. Jiittner 
derived the thermal equilibrium distribution for a rela- 
tivistic gas already in 1911 fl3 |. After that many au- 
thors have given important contributions to the develop- 
ment of the relativistic kinetic theory (for an introduc- 
tion see e.g. [l5j],[l(|). In spite of this progress in recent 
years a controversial debate about the correct general- 
ization of Maxwells velocity distribution in special rel- 
ativity arose [H, [T3|- 0- Recently numerical micro- 
scopic l-dimension"simuktions and a critical anal- 
ysis of alternative findings [2~fl | yield arguments in favor 
the Jiittner distribution. Relativistic diffusion processes 
has been comprehensively reviewed 11| which besides 



the issue of stochastic relativistic diffusion theory also in- 
cludes relativistic equilibrium thermostatics and micro- 
scopic models for Langevin-typ equations and where a 
more complete list of references can be found. 

Diffusion theory in the Euclidian space R d is a well 
developed topic (see e.g. (22[, [3]). However, the de- 
scription of diffusion on a non-Euclidian manifolds M d is 
a subject containing several pitfalls. There exists a well 
developed rigorous mathematical theory of stochastic dif- 
ferential equations and diffusion processes on Rieman- 
nian manifolds with a definite metric signatur (see e.g. 
[IS] , [IH ) • The stochastic calculus on Riemannian mani- 



folds found considerable interest in mathematics and has 
played a central role in recent years within the analysis in 
path and loop spaces in topology and other fields. How- 
ever, this mathematical approach can not be applied to 
describe diffusion on Lorentzian or Pseudo-Riemannian 
manifolds with an indefinite metric. In the present pa- 
per we derive a physically motivated modification of this 
calculus to describe diffusion in the phase space of the 
position and velocity, in which the difficulties in the de- 
scription of diffusion on manifolds with indefinite met- 
ric signature are bypassed. In this approach a general- 
ized relativistic Langevin equation in the fiber bundle of 
position, velocity and orthonormal velocity frames is de- 
fined from which the generalized non-relativistic Kramers 
equation in external force fields is derived. As will be 
shown the derived relativistic diffusion equation satisfies 
the general principle of special relativity and is invariant 
under Lorentz transformations. The steady-state solu- 
tion of this equation for a heath bath with constant fric- 
tion coefficient yields the Jiittner distribution. 
The paper is organized as follows. In Chapter 2 the con- 
cept and main formulas of the mathematical stochastic 
calculus on Riemannian manifolds is presented. Since 
this approach is little-known in physics in the appendix 
B some details of this calculus and its relation with the 
stochastic calculus on Euclidian manifolds (appendix A) 
are presented. In Chapter 3 the generalized relativistic 
Langevin equation in the fiber bundle of position, veloc- 
ity and orthonormal velocity frames is defined and the 
relativistic diffusion equation for the probability density 
function or the transition probability is derived. In Chap- 
ter 4 the steady-state solution for a relativistic gas in a 
heat bath with constant friction coefficient and in Chap- 
ter 5 the non-steady solution for the force-free case are 
derived and in chapter 6 the conclusions are presented. 
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2. MATHEMATICAL STOCHASTIC CALCULUS 
ON RIEMANNIAN MANIFOLDS 

Stochastic differential equations in diffusion theory 
in a d-dimensional Euclidian space R d with continuous 
pathway are defined by the fundamental d-dimensional 
Wiener process W a (t). On a Riemannian manifold M d 
the fundamental Wiener process is difficult to handle. 
By using an inadequately posed formulation of a stochas- 
tic differntial equation it is not assured that its solution 
remains on the manifold M d which leads to inconsis- 
tent results. The key idea in the mathematical concept 
of diffusion on general d-dimensional Riemannian man- 
ifolds M d (with definite metric signature) is to define a 
stochastic process on the curved manifold using the fun- 
damental Wiener process each component of which is a 
process in the Euclidian space R d [2ij],[2^]. Intuitively, 
we can understood this concept as follows. Consider 
a two dimensional stochastic motion of a particle on a 
plane. If the trajectory of the particle is traced in ink 
and a sphere on the plane is rolled along the stochas- 
tic curve without slipping the resulting path which is 
thus transferred defines a random curve or a stochas- 
tic Markovian process on the sphere. This method can 
be applied for diffusion on a general Riemannian mani- 
fold. The tangent space of a Riemannian manifold is en- 
dowed with an Euclidian structure and therefore we can 
move the manifold in the tangent space by construction 
of a parallel translation along the stochastic curve with 
the help of the orthonormal frame vectors e a = e l a (x)di 
(i,a=l...d) and the Christoffel connection coefficients 1^, , 
x = [x\...Xd)i ®i — djdx 1 . In the common mathematical 
denotation this means that the random curve on the Rie- 
mannian manifold M d is lifted to the horizontal curve on 
the frame bundle O(M) and this horizontal curve corre- 
spond uniquely to a random curve in an Euclidian space. 
In local coordinates on a Riemannian manifold the in- 
finitesimal motion of a smooth curve c z (t) in M d is that 
of 7* (t) in the tangent space (which can be identified with 
R d ) by using a parallel transformation: dc l = e l a (x)d"f a 
and de l a {x) — —F l ml e l a dc m . Therefore a random curve 
can be defined in the same way by using the canonical 
realization of a d-dimensional Wiener process (defined 
in the Euclidian space) and substituting d-f a — > dW a {t). 
As explained in the appendix B the stochastic differential 
equation describing diffusion on a Riemannian manifold 
in the orthonormal frame bundle O(M) with coordinates 
O(M) ={x l ,e l a } then is given by 

dx l {r) = e^rjodW? + A l (r)dT, (1) 
de^r) = -r ml e l a odx m (r) 

Here r ? o6 et(x(r))4(x(r)) = fl « ftej = -I^e*, 

is the Riemannian metric and i] ab — 5 a the flat 

Euclidian metric where S ab is the Kronnecker sym- 



bol. The components of the elementary Wiener pro- 
cess dW a = W a (t+At)-W a (t) are defined in the Eu- 
clidian space with the probability density P(W a ) = 

(2DirAt)~^ exp(— ^jj^t ) an< ^ with the expectation 
values (W a ) = and (V°(t)W 6 (t + s)) = DsS ab . 
A 1 are the components of an arbitrary tangential vector 
and D is the diffusion constant which here is independent 
on the time and space variables. Eq. (1) is defined in the 
Stratonovich calculus (denoted by the symbol o). 

Associated to each diffusion process, there is a second 
order differential operator denoted as the generator A of 
the diffusion. This operator is associated with the Kol- 
mogorov backward equation and defined in appendix A 
and B for diffusion processes on Euclidian and Rieman- 
nian manifolds, respectively. For the stochastic process 
lifted to the fiber bundle O(M) as defined in Eq.(l) the 
diffusion generator in the Stratonovich integral interpre- 
tation is given by Eq.(B4) as 



a O(M) = j H a H a + H (2) 

a=l 

where the fundamental vector fields H a and Hq are given 
by Eq.(B3): 

H = A i {r,X)d i -T i ml e l a A m {r)^ 

The diffusion generator A (m) on O(M) can be pro- 
jected to M d with f(r) = /(x,0), r = (ir l ,e*) using 
the relation A G ( M )/(r) = Am/(x) where 

3 

Am = f 5>* fite&O + A% = (—Am + A%) (4) 

a=l 

and Am = g^didj — g^T^dk is the Laplace-Beltrami 
operator on the manifold M d . The generalized Fokker- 
Planck equation is obtained by the adjoint of the diffusion 
generator A^ (which includes the volume element ^/g, 
g = det{g iJ }). Since the Laplace-Beltrami operator is 
self-adjoint Am = A^ f this equation takes the form: 

<9$ D 

— = -div x (A$) + -A M $, (5) 

UT Z 

where div x (A&) = g~^ di(g^A l $) is the divergence op- 
erator on the Riemannian manifold, <f> = $(x, r | y,0) 
is the transition probability with the initial condition 
<&(x, | y, 0) = S(x-y) and adequate boundary condi- 
tions at infinity. The probability density Lp(x l ,r) is de- 
termined by the same equation with the initial condition 
p(x,r = 0)=¥>°(x). 
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A remarkable particularity of Markovian diffusion on a 
Riemannian manifold is the supposition that for the dif- 
fusion coefficients in Eq.(l) only the orthonormal frame 
coefficients e„(x) are admissible which are directly re- 
lated with the geometry of the Riemannian manifold. In 
contrast on an Euclidian manifold a much more general 
class of diffusion coefficients is permitted. 



3. RELATIVISTIC DIFFUSION IN THE PHASE 
SPACE 

A direct application of the mathematical calculus of 
diffusion processes on general Riemannian manifolds for 
relativistic physics is not possible due to the supposition 
restricting the stochastic formalism to the special case 
of a Riemannian manifold with a definite metric signa- 
ture, but in relativity theory the Lorentzian or Pseudo- 
Riemannian manifolds exhibit an indefinite metric sig- 
nature (-,+,+,+). Moreover, it has been proven that 
Markovian diffusion processes in the base manifold (posi- 
tion space) on a Lorentzian or Pseudo- Riemannian mani- 
fold do not exist However, considering more care- 
fully the mathematical model described by Eq.(l) one 
can recognize a significant difference to the physical non- 
relativistic diffusion model described by the Langevin 
equation. In Eq.(I) the noise term directly acts on the 
position variable in the base space, while the noise term 
in the Langevin equation operates like a force on the 
change of the velocity in the tangent space. This crucial 
difference in the mathematical model to the physically 
motivated Langevin approach is the central point which 
enables a generalization of the Markovian diffusion the- 
ory within the framework of the special relativity theory 
performed in the phase space of coordinates x 1 and nor- 
malized velocity variables u % (i = 1,2,3). The 4-velocity 
space for a gas with massive particles is a hyperboloid (or 
pseudo-sphere) described by the relation (u ) 2 — (u 1 ) 2 — 
(u 2 ) 2 — (u 3 ) 2 — I. This means that the relativistic veloc- 
ity space is a noncompact hyperbolic 3-dimensional Rie- 
mannian manifold (and not Pseudo-Riemannian). There- 
fore with an appropriate modification we can apply for 
the velocity space the mathematical stochastic calculus 
for Riemannian manifolds as presented in the appendix 
B. Since the stochastic force acts directly only on the 
change of the velocity and not on the position coordi- 
nates we can define the relativistic generalization of the 
Langevin equation in a fiber bundle here denoted by 
F(M L ) :{x\ u\E l a } = F{M L ) where x i belong to the 
Lorentzian base manifold Ml, the relativistic velocity u l 
to the tangent space TMl and E % a (u) are the moving or- 
thonormal frames in the hyperbolic velocity space (this 
means they belong to the second order tangent space 
TTMl). Locally, this fiber bundle is simply the product 
space of these three sub-spaces. With a corresponding 
modification of Eq. (Bl) in appendix B the generalized 



relativistic Langevin equations can be defined in the fiber 
bundle space F(Ml) by 

dx l (r) = M*(r)dr, 

du l {r) = El{T))dW a {T) + F l {T)dT, (6) 
dEl{r) = -f ml (u)E l a {T)odu m (r) 

= -Y ml ^)E l a (r)[Er(r))dW b (r) + F m (r)dr}. 

Here r is an evolution parameter along the world lines of 
the particles which can be chosen as the proper time. 
The laboratory time t — tu°/c is a function of the 
proper time r and u° which here and below is defined 
by u° = [1 + {u 1 ) 2 + (u 2 ) 2 + (u 3 ) 2 ]i . 7* n; (u) are the 
Christoffel connection coefficients on the hyperboloid and 
F l =K l /m, where K l are the spatial components of the 
4-force, m is the rest mass of the particles and the indices 
a,b denotes the spatial components in the hyperbolic 
space (a,b=I,2,3). Since the stochastic force dW a (r) do 
not act directly on the position variable x 1 (t) the in- 
definite signature of the Lorentzian manifold here do not 
create any difficulty, as it arise for a stochastic differential 
equation as Eq.(l) for a manifold with indefinite metric. 
Sufficient conditions for the existence and uniqueness of 
the stochastic differential Eq.(6) are that the drift and 
diffusion coefficients satisfy the uniform Lipschitz con- 
dition and the stochastic process X(t)={x(t), u(t)} is 
adapted to the Wiener process W a (r) , that is, the output 
X(r2) is a function of W a (ri) up to that time (r% < t<z) 
[23j . The moving frames in the hyperbolic velocity space 
are defined by the relation 

a=l 

or equivalently 

G l3 KEi=5 ab , (8) 

where Gy is the Riemannian metric of the hyperbolic 
velocity space, G u the inverse matrix of Gy and the 
Christoffel connection coefficients 7„/(u) on the hyper- 
boloid are given by 

7}fe(u) = ^G im [dG jm /du k + dGmk/dv? - dG jk /du m }. 

(9) 

Since the manifold on the hyperboloid is embedded into 
the Minkowski space, the metric Gij(u) can be calcu- 
lated from the infinitesimal arc length given by ds 2 = 
-(d M °) 2 + (dM 1 ) 2 + (du 2 ) 2 + (du 3 ) 2 with u° = [l + (u 1 ) 2 + 
(m 2 ) 2 +(m 3 ) 2 ] 2 . In this way we obtain ds 2 = Gy {i^difdui 
with Gij(u) = Sij - {uiUj)/^ ) 2 ^ = detGy = (w )~ 2 
and 7j fe (u) = —u l Gjk- Corresponding the definition of 
fundamental vector fields on 0{M) in the appendix B 
one can now introduce the fundamental horizontal vec- 
tor field H a and H on the fiber bundle F(M L ). With 
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corresponding modifications we find from Eq.(B3): 

d_ 



d 



(10) 



and the diffusion operator Ap(M L ) f° r the stochastic 
process is given as in Eq. (B4) by 



In the relativistic framework, the Lorentz-invariance of 
the physical laws is one of the most fundamental prop- 
erty. Eq.(13) refers to a special inertial rest frame X of 
an observer. Now consider a second observer at rest in an 
another inertial frame X that moves with constant veloc- 
ity w relative to X. As shown and discussed by many au- 
thors (see e.g. ll|,[l||) the probability density function 
<P(t; x, u) transforms as a Lorentz scalar; i.e. it fulfills 
the condition 



= — 53 HaHa 



L F(M L 



(11) 



We project the stochastic curve from the fiber space 
F(Ml) with coordinates r ={x l , u % , E l a } to the phase 



space with coordinates {x\u 1 } 



^F(M L ) 



m = 



Ap/(x, u, 0) , where the diffusion generator in the phase 
space Ap is given by 



D ^ i 



+ u l d/dx l + F l d/du\ (12) 



The special relativistic diffusion equation in the phase 
space is given by the adjoint of the operator Ap and 
similar to Eq.(B7) the generalized relativistic Kramers 
equation takes the form 



8$ 9$ d 

T*"*' (13) 

where A u is the Laplace Beltrami Operator of the hyper- 
bolic velocity space given by 



A„ = G 13 



if- 



dvtdui 



(14) 



and the corresponding divergence operator is given by 



(15) 



with G = det {G i3 }. 

Eq. (13) represents the relativistic generalization of 
the Kramers equation for the probability density func- 
tion $ = 4>(t; x, u) with the initial condition 4>(t — 
0;x, u) = </> (x, u). The transition probability is de- 
termined by the same equation, but is defined by the 
initial condition $(x, u, | xo,uo,0) = ~^^( ul ~ u h) 
5(u 2 - u 2 )S(u 3 - ul)S(xi ~ x° 1 )S(x 2 - x° 2 )S(x 3 -x° 3 ).If the 
force F l depend on the time t we have to substitute t by 
t = r[l + (u 1 ) 2 + (u 2 ) 2 + (u 3 ) 2 ]3/c. For an external elec- 
tromagnetic field F^ v the normalized force F l is given 
by F l = eFlu v . For small velocities (| u l | 2 < 1) Eqil3) 
pass over to the nonrelativistic Kramers equation [26j |. 



0'(r',x',u') =0(r,x,u) (16) 
where the variables r',x', u' 

x" = A)x j + A^x°, u' j = K)u 3 + A*u°, r' = r (17) 

are related by the Lorentz transformation. This re- 
quires that Eq.(13) is invariant with respect to a Lorentz 
transformation. By using x — tm°, the chain rule 
d/dx z = K^d/dx 3 and the inverse transformation u 1 — 

A)u 13 + A* u'with A^A^ = 6% we find 



U °A* 



dx l 



dx l 



dx l 



(18) 



df 



dx l 



On a Riemannian manifold the divergence operator 
and the Laplace-Beltrami operator are intrinsically in- 
variant with respect to general coordinate transforma- 
tions. Therefore Lorentz transformation on the pseudo- 
sphere do not change the explicit form of these oper- 
ators. For the divergence operator this can be simply 
proven by the transformation property of the covari- 
ant differentiation DjF 1 = 8F l /duj + ~l l j k F k given by 
D'jF li = (du li /du k )(du l /du l3 )D l F k . For the diver- 
gence operator div u (F(j)) = Dj(F 3 (j)) this yields the re- 
lation D' jF'J =DjF 3 . The invariance of the Laplace- 
Beltrami operator A u <j) = Dj{G^ di4>) can be similarly 
proven. Since G v di4> = A 3 transforms like a vector and 
the divergence DjA J is as shown above an invariant oper- 
ation the relation A u (f> = A u <cj) follows. By using Eq.(18) 
we express the variables x l , u 3 by the new variables 
x n ,u 13 in the inertial frame X and account the invari- 
ance of the divergence and Laplace-Beltrami operator. 
Then Eq.(13) takes the form 



iJ^L-div^F'tf + jAv't, (19) 



Thus, the derived relativistic diffusion equation satis- 
fies the general principle of special relativity and is invari- 
ant under Lorentz transformations. Note that Eq.(13) 
differs from previously derived relativistic diffusion equa- 
tions. Debbasch et al. 6] introduced a phenomenological 
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relativistic Langevin equation in the phase space and de- 
rived from this a generalized Kramers equation of the 
classical Ornstein-Uhlenbeck process in which the diffu- 
sion term is given by the 3-dimensional Euklidian Lapla- 
cian in the momentum space. A different approach has 
been presented by Dunkel and Hanggi § but the final 
generalized Kramers equation and its steady-state solu- 
tion also differ from Eq.(13). Both the derived equa- 
tion in [6] as well those in in Hare not invariant under 
Lorentz transformations (see [ll| and [27]]). 

For the solution of the relativistic diffusion equation 
it is convenient to introduce the hyperbolic coordinate 
system for the 4- velocity defined by u 1 = shasin^sint/j, 
u 2 = sh a sin z? cos tp, u 3 = sha cos $ and u° = cha. We 
denote the velocities in the non-Cartesian coordinates 
by u 1 = a,u 2 — 9,u 3 = p ,a = 1,2,3. The metric in 
this coordinates are simply to calculate and are given by 
Gn = 1, G22 = sh 2 a,Gs3 — sh 2 a sin 2 tp and Gij = 
for i 7^ j. With the given metric the Laplace Beltrami 
Operator A„ in the hyperbolic velocity space takes the 
form 



A = 



da 2 



(sha) 



2ctha- 

oa 

d 2 „ d 



(20) 



& 2 



(sin i?) 2 dip 2 



and 



div u (F$) = (sh 



(21) 



2 jL(( sha ) 2 F «<I> 
oa 

-(shay^sintf)- 1 -^- (smi?F"$) 
ov 

-(sha)- 1 (sin^)- 1 ^-( J F*'$), 
dip 

is the divergence operator in the hyperbolic veloc- 
ity space. Here the force in the hyperbolic coordi- 
nate system F a , F® , F v is related with F i by F a = 
(cha)" 1 [sin^cos^ 1 + sin^F 2 ) + costfF 3 ], F" 6 = 
(sha) -1 [cosi^cos^-F 1 + sinpF 2 ) — sin^F 3 ] and F v — 
(sha) -1 (sint9) -1 [— s'mipF 1 + cos pF 2 ] with the ini- 
tial condition 0(r = 0; Xi, a, tp) = 4> (xi, a, ip). 
The transition probability is determined by the ini- 
tial condition ^(x, a, 1?, tp, j xo, ao, i?0) Vo' 0) = 
(sha) -2 (sini9) -1 <S(a - a ) S(d - & )6(tp - p )S( Xl - 
x 1 )S{x2-x°)5(x 3 -x° 3 ). 

The diffusion of massless particles such as e.g. the 
diffusion of photons in random media (see e.g.[30j]) can 
be described in analogous way, but with the condition 
{u ) 2 - (u 1 ) 2 - (u 2 ) 2 - (u 3 ) 2 = 0. For the velocity co- 
ordinates now we choose u° — r, u 1 — r sin^sintp, u 2 — 
r sin ■& cos tp, u 3 = r cos $. The Laplace-Beltrami operator 
then takes the form 

d 2 d 1 / d 2 d Id 2 

-2r- - — * 



A = 



dr 2 



dr 



a d 



(sini?) 2 dp 2 

(22) 



4. STEADY STATE SOLUTION OF PARTICLES 
IN A HEAT BATH: THE JUTTNER 
DISTRIBUTION 

First, we consider particles with a rest mass m in a 
gas in an isotropic homogenous heat bath. The interac- 
tion of particles with the bath is described by a random 
noise force and a friction force. In the nonrelativistic 
theory the friction force is given by /* = -cmv 1 where 
v is the friction coefficient and v 1 are the components 
of the nonrelativistic velocity. The relativistic general- 
ization of the friction force requires the introduction of a 
friction tensor v % a similar to the pressure tensor in special 
relativity theory Q, Q. The friction force is expressed 
as F l — v l a [u a — U a ], where U a is the 4- velocity of the 
heat bath. For an isotropic homogeneous heat bath the 
friction tensor is given by 



< = v{ri i l3 + u i u a ), 



(23) 



with v denoting the scalar friction coefficient measured in 
the rest frame of the particles. In the laboratory frame 
the heat bath is at rest described by U a — (1,0,0,0). 
Therefore the friction force is given by F l = —vu l vP or in 
hyperbolic coordinates F a = — vsha,F® = F v = 0. We 
consider the spatial homogenous and isotropic solution 
of equation Eq.(13) with Eq.(20) and (21) described by 



dc/)(a) 
dr 



= D 



d 2 

da 2 



2— etna 
oa 

d 



(«) + 



v(sha)- 2 — r(sha) 3 0(a)l . (24) 
Or* 



The steady-state solution of this equation is given 
by (f>(a) — Cexp{— xcha} with \ — v/D and 
C =4irK2(x)/x- ^(x) denotes the modified Hankcl 
function. This distribution is identical with the Juttner 
equilibrium distribution if the Einstein relation kT = 
Dmc 2 /v and consequently % = is used. The above 
derived relativistic diffusion equation Eq.(13) yields for 
the 3D case the correct thermodynamic relativistic equi- 
librium distribution for a constant friction coefficient. 
Note that in a previously derived relativistic diffusion 
equation [f| the Juttner equilibrium distribution for a 
relativist gas only arise as the steady-state solution 
for a specifically adapted energy-dependent friction con- 
stant v = vq (u ) 2 . The generalized Kramers equation 
in [H yields for constant friction coefficients and in the 
Stratonovich interpretation a modified Juttner distribu- 
tion, but using the non-standard pre-point discretization 
rule the Juttner distribution was found. Recently fully 
relativistic one-dimensional molecular dynamics simula- 
tions favored the Juttner distribution in the ID case [Til ] . 
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5. NONSTEADY SOLUTION FOR THE 
FORCE-FREE CASE 



Now we consider the unsteady solution of Eq.(13) 
for a spatial homogenous gas with vanishing force F\ 
We use the Laplace transformation 1?, ip, r) = 

oc 

/ $(a, d, <p, A) exp(— Xr)dX and for the eigenvalue func- 
o 

tions the ansatz <& J M (a, ip, A) = gj(o;)Y^(&, y>) where 
Yjtf{0, f) —V J M (d)e lMv are the spherical harmonics with 
the associated Legendre functions ¥ J M (■#). For <?}(«) the 
following eigenvalue equation is derived: 



D 



d 2 
da 2 



d , J J+ 1 

2— ctha ^5 - 

eta sh a 



\\g$(a)=0. 



(25) 

Here the discrete index J takes the values J = 0, 1, 2, . . . 
and M = — J, — J + 1, ... 0, 1, ... J. The eigenfunctions 
for the Laplace-Beltrami operator with the eigenvalues 
A = D(x 2 + 1) satisfying the boundary condition are 
given by 



gj(z) = 2 - 1) 



J+i 



cos(x archz), (26) 



with z = cha and = (-1)' 7+1 fl -/M^ 2 + k 



fe=0 

The eigenfunctions $}[ / (a, y>, x) satisfy the relations of 
orthogonality and completeness. The transition probabil- 
ity is determined by the initial condition $, ip, t — 
| a o ,0 o ,cp o ,O) = (sha) _2 (siiii?)- x <$(a - a ) 5(0- 
$0)8 (ip — ip ). Using the orthogonality relation we can 
write 



3>(a, 0, ip, t 



a ,do,P ,0) (27) 
$ J M (a, 0, V , x)3>*/ (a , a ,<p o , x) 



E 

M,J 



cxp[-L>(x 2 + l)r]dx. 

Let us now consider the fundamental solution J = 0. 
Substituting 3^ (a) = — -^-(sh a) -1 s'mxa into Eq.(27) 
gives the transition probability 



$(a, t I a , 0) = C(tD) _ 3 exp{-r.D} 



x exp 



sh a sh «o 

1 (28) 



a 2 + a 2 
4Dt 



substituting in Eq.(25) sha — > a,cha — > 1 we find for 
J=0 go (a) = — -^=^(a) _1 sinxa, but the eigenvalue is 
given by A = Dx 2 and the solution now is 



*(a, t I a , 0) = C(Dr)~^a- 1 ao 1 



exp 



exp 



(a - ap) 2 
4L>t 



(a + ap) 2 



(29) 



For «o = Eq.(29) pass to the known Wiener distri- 
bution in the velocity space. In the limit a -C 1 the 
short time behavior of Eq.(29) is up to an exponential 
small factor in agreement with Eq.(28), however in the 
long time behavior both solutions differ by the exponen- 
tial factor exp(— Dt). In Fig.l the relativistic distribution 
Eq.(28) is presented by the solid lines for ccq — and the 
Wiener distribution Eq.(29) by the dotted lines. As can 
be seen both distributions differs by orders of magnitudes 
even in the non relativistic region »« 1 for long times 
Dt 3> 1. This discrepancy can be explained by the topo- 
logical properties of the hyperbolic space (included by 
the boundary conditions) which are different from that 
of the Euclidian space in the non relativistic theory. The 
deep connection between local and global properties of 
diffusion processes is a central topic in the mathemati- 
cal field of heat kernels on Riemannian manifolds [29| . 
The appearance of the factor exp(— Dt) can also be ex- 
plained by physical arguments, it comes from that in the 
hyperbolic coordinates the Jacobian is proportional to 
sh a which is exponentially large for r — > 00. For large r 
the entire velocity space is explored and the small factor 
exp(— Dt) cancels the exponentially large Jacobian and 
guarantees the probability conservation. 

Let us still compare the solution for a massive par- 
ticle with that of a massless one. With the Laplacian 
Eq.(22) the eigenvalues takes the form A = Dx 2 and the 
eigenvalue solutions are $ J M (a, 0, ip, A) = gj(a)Y^(6, (p), 
gj(z) = ^(^(-ly^d-Y^). The photon tran- 
sition probability for J = then is given by 

*(o,t I oo,0) = (2rr V^)- 1 (rD)-hh(^-) 

ZJJt 

ex P {-^}. (30) 

In comparison with Eq.28 the factor exp(— Dt) here is 
absent since the Jacobian do not increase exponentially 
for r — > 00. 



6. CONCLUSIONS 



with C = 2(47r)~5. For ao — ► this solution was first 
found in [28j]. For small velocities (a « 1) the transition 
distribution shows a remarkable behavior. If we solve 
the corresponding non-relativistic Kramers equation [26| 



In conclusion, a theory of Markovian diffusion pro- 
cesses within the framework of the special theory of rel- 
ativity is formulated. In the derivation of the basic rela- 
tivistic diffusion equation the mathematical calculus of 
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a 



FIG. 1: Relativistic (solid thick lines, Eq.(28)) and non rela- 
tivists (dotted thin lines, Eq.(29)) distributions for different 
times Dt and ao — 

stochastic differential equations on Riemannian mani- 
folds is used, which here is modified for the description of 
diffusion in the phase space on Lorentzian manifolds with 
an indefinite metric. A generalized Langevin equation in 
the fiber space of position, velocity and orthonormal ve- 
locity frames is defined and the generalized relativistic 
Kramers equation is derived. This diffusion equation is 
invariant under Lorentz transformations. In the case of a 
relativistic gas in a heat bath its steady-state solution is 
identical with the Jiittner distribution for constant fric- 
tion coefficients. An analytical nonsteady solution for 
the transition probability is given for the special case of 
vanishing external fields. This solution differs from the 
Wiener velocity distribution even for small velocities due 
to topological reasons. 

The formalism presented in this paper can be extended 
to a theory within the framework of general relativity. 
This will be done in a forthcoming paper. 

APPENDIX A: DIFFUSION ON EUCLIDIAN 
MANIFOLDS 

Using appropriate modifications the stochastic calculus 
on Riemannian manifolds utilizes basic theorems and for- 
mulas from the calculus on Euclidian manifolds. There- 
fore in appendix A some basic notations, formulas and 
theorems for the Euclidian space are summarized includ- 
ing some elementary proofs of the basic theorems. 

Diffusion processes in a d-dimensional Euclidian space 
are described by stochastic differential equations of the 
form [H- [H 

dX l = a l a {T,X)dW a + b i (T,X)dr. (Af) 
X = (Xi,...,Xd) G R d is a stochastic process with 



X(0) = x; x = (x^...x<i) with the time r (t > 0). The 
diffusion coefficients <t^(t, X) are given matrices and 
the drift coefficients 6*(r,X) are coefficients of a smooth 
vector field. W a are the components of the elementary 
Wiener process dW a = W a (t+At)-W a (t) with the proba- 
bility density P(W Q ) = {2itDAt)-? exp(-^g^) and 
the expectation values (W a ) = 0, (VK a (r)WV + s)) = 
DsS a t,, D is the diffusion constant. Eq.(Al) can be trans- 
formed into an integral equation 

T T 

X l r ^X^ + J a l a (s,X)dW? + J b l (s,X)ds. (A2) 
o o 

The stochastic integral in the second term of 

T 

Eq.(A2) is defined as the limit f <r t a (s,X)dW? = 

o 

n 

E o-i(s*,X)(W a (si) - W a (si_i)) as n -► oo. This in- 
»=i 

tegral depends on the choice of the intermediate point 
s*. With the choice s* — s.j_i (post point rule) the Ito 
stochastic integral is defined. The Ito integral is a Marko- 
vian process and plays a fundamental role in the theory 
of diffusion processes and most of general mathemati- 
cal treatments can only rigorously proven by using this 
calculus. Alternatively, choosing s* — Sj_i (mid point 
rule) the Stratonovich stochastic integral is defined. The 
Stratonovich integral has the advantage of leading to or- 
dinary chain rule formulas under a transformation. This 
property makes the Stratonovich integral natural to use 
for stochastic differential equations on Riemann mani- 
folds. However, in general Stratonovich integrals are not 
Markovian processes which hinders rigorous mathemati- 
cal treatments in most cases. 

A differential equation Eq.(Al) is defined only with 
respect to one of the both stochastic integrals, changing 
the interpretation of the integral refers to a differential 
problem with different solutions. Therefore the chosen 
interpretation should be denoted in the differential equa- 
tion. The symbol cr l a {r, X)dW a denotes the Ito integral 
interpretation and cr*(r,X) o dW£ the Stratonovich in- 
terpretation. 

With the Ito interpretation the solution X l T of Eq.(Al) 
is denoted as an Ito process if the diffusion and drift co- 
efficients satisfy the Lipshift condition and cr^Cr-X) is 
adapted to the fundamental Wiener process [23j . An 
Ito process has the important property of being Marko- 
vian. Then Y T = f(X T ) is also an Ito process, and its 
stochastic differential equation (Ito formula) is 

d f = §£k( t > x k( t > x )w+ 

a=l 

+ b l (r, X)d.J}dT + Da l a (r, X)d l fdW a . (A3) 
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Associated to an Ito process is the diffusion generator 
A of X r , which is defined to act on a suitable function / 
by: 



A f = lim 

t-»o 



£*[/(X T )] - /(x) 



t 



(A4) 



where x=X Q is the initial point of X T . By using Eq.(A3) 
one can show that A is given by 



f d 

/ f(*)-Q^4 l (T,X-)dx 1 ..dx d 



(A10) 



We integrate by parts and discard surface terms to obtain 



J f(x)A*(j)dxi..dxd = 



A f= f E4(T J X)ai(r J X)^/ + 6 i (r J X)ft/. (A5) 

a=l 

The generator A can be used in the derivation of Kol- 
mogorov's backward equation. This equation describes 
how the expected value E x [f(X T )] of any smooth func- 
tion / of X evolve in time. If we define 



u(r,x) = B»[/(X T )] 
then u satisfy the following equation: 



(A6) 



f(x)—(/)dxi..dx ( i, 



(All) 



where A*denotes the Hermitian adjoint of A. Since /(x) 
is arbitrary we find for the Fokker-Plack equation within 
the Ito integral interpretation: 



dr 



(r,x)=^*0(r,x) 



with the adjoint operator 



(A12) 



9 u(r,x) = E*[±f(X T )] = ^E*[f(X r )] 



dr ' dr 
and with Eq.(A3) 



d_ 



u(r, x) = Au(r, x) 



(A7) 



(A8) 



with w(0,u)=/(x). The Fokker-Plack equation (or for- 
ward Kolmogorov equation) describes how the proba- 
bility density function </>{t, x) of X r evolve with time. 
This density function is defined as a function of the d 
stochastic variables, such that for any domain Q in the 
d-dimcnsional space of the variables X\, Xd, the prob- 
ability that a realization of a set of variables falls inside 
the domain CI is 



Pr(A%\...A^ en) = J </>(T,x)dx 1 ...dx d . 



(A9) 



The probability density function can be used to cal- 
culate the expected value E x [f(X T )} by E x [f{X T )} 

=J f(x.)<fi(r,x)dxi...dxd- We now consider the time 
Q 

derivative of such expectation value and use again Ito's 
formula Eq.(A3): 



A*f=j^d i d j ai(r,X)ai(T,X)f-d i b i (T,X)f. 

a=l 

(A13) 

Since the stochastic calculus on Riemannian manifolds 
is naturally formulated in the Stratonovich integral inter- 
pretation, we will consider the connection between both 
types of integrals. Let us formulate the stochastic dif- 
ferential equation Eq.(Al) with the Ito interpretation by 
an corresponding equation with the Stratonovich inter- 
pretation 



dX l = 5* (r, X) o dW a + 6*(r, X)dr. 



(A14) 



In the Stratonovich interpretation the differentiation of a 
function / yields the ordinary chain rule, e.g. in the Ito 
formula Eq.(A3) the first term with second-order deriva- 
tives does not appear and the above given treatment can 
not be done in a consistent way. However, there exist 
a connection between the Ito and the Stratonovich inte- 
grals [22 



25j of functions y> a (X(r),r) in which X(r) is 



the solution of the Ito differential equation Eq.(Al): 



T T 

J Va { S ,x)dw: = J <p a ( Sl x)odw:- 



E x [±-f(X T )] = j-E x [f(X T )) = E*[Af(X T ) 



J <p(r, x)Af(x)dx 1 ..dx d 



D 
~2 



T 

[ ai(T,X)dMs,X)odW° (A15) 



If we now make the choice 
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d<(x) 



■ ml a 



dx r ' 



(Bl) 



#(t,X) = 6<(r,X) - ^ ^ (r, X)^<(r, X), 



<(r,X)=^(r,X) 



(A16) 



and substitute 6 l (r, X) into Eq.(A5), the diffusion oper- 
ator A in the Stratonovich interpretation is 



A = f £<(r,x)9 i( ri(T,x)^ + b i {r,x)d i . (A17) 

a=l 

The Fokker-Plack equation Eq.(Al2) with respect to the 
Stratonovich interpretation is then given by: 



d D . d . , 

a— 1 



-^(t.x^t.x). (A18) 
We Introduce the fundamental vector fields 



L a = a* (r,x)9i, L = b l (T,x)d i 



(AI9) 



the generator A of the stochastic process can be ex- 
pressed by the fundamental vector fields L a 



Here ?7 V (x)e l a (x) = ^,c^ a = -I*e», ff y is the Rie- 
mannian metric and ?7 ab = <5 ab the flat Euclidian met- 
ric where S ab is the Kronecker symbol. r l (t) is called 
the horizontal lift of the curve x % (t) to the orthonormal 
frame bundle O(M) and it lies in the Euclidian space 
R d+d _ The horizontal 

curve 7 l (i) corresponds uniquely 
to a smooth curve in the tangent space (which can be 
identified with an Euclidean space R d ). 

Stochastic differential equations on a Riemannian 
manifold can be defined by using the above described 
horizontal lift to the orthonormal frame bundle 0{M) 
endowed with an Euclidian structure 24j , 25| . By using 
this approach a stochastic process on the Riemannian 
manifold can be constructed by using the fundamental 
Wiener process each component of which is a process 
in the Euclidian space R d and interpreting the corre- 
sponding stochastic integral in the sense of Stratonovich. 
This means that the manifold is moved along a stochas- 
tic curve in the tangent space by a parallel translation 
with the help of the orthonormal frame bundles and 
the Christoffel connection coefficients Y 3 ib . Correspond- 
ingly a random curve can be defined in the same way 
as in Eq.(BI) by using the canonical realization of a 
d-dimensional Wiener process and substituting d-y 01 — > 
dW a (t). Therefore the stochastic differential equation 
describing diffusion on a Riemannian manifold is 

M 



dx i = ej,(r)o dW a + A i dr, 



a=l 



L . 



(A20) 



APPENDIX B: DIFFUSION ON RIEMANNIAN 
MANIFOLDS 

Stochastic differential equations are defined by the 
driving Wiener process W a ( or more general by a semi- 
martingale), but this process is difficult to handle on a 
Riemannian manifold. In differential geometry for a gen- 
eral d-dimensional Riemannian manifold M d (with def- 
inite metric signature) equipped with a Christoffel con- 
nection Y J ib it is possible to lift a smooth curve c 1 (t) in M d 
to a horizontal curve in the tangent bundle TM which is 
endowed with an Euclidian structure by using the bun- 
dles of orthonormal frames e a — e l a (x)di (i,a=l-d). The 
orthonormal frame bundle 0(M ) is described by the lo- 
cal coordinates {r = (x l , e*)} = 0(M). The infinitesimal 
motion of a smooth curve x l (i) in M d is that of j l (t) in 
0(M) described by the ordinary differential equations for 
a parallel transport 

dx* = <(x)d 7 Q , 



dei(r) = -r ml e l a (r)o^ 



(B2) 



where the components of an arbitrary tangential vec- 
tor A 1 are additionally introduced for a more gen- 
eral situation with account of an external force field. 
The components of the elementary Wiener process 
dW a = W a (t+At)-W a (t) are defined in the Euclid- 
ian space with the probability density P(W a ) = 

(2DAt)~i exp(— ^J^t ) w * tn tne expectation values 
(W a ) = 0, (W a {T)W b (T + s)) = Ds5 ab . 

The derivation of the Kolmogorov backward equation 
with the definition of the diffusion operator Aq(m) can 
be performed by the same rules as in Euclidian space in 
the Stratonovich calculus. Corresponding the definition 
of the fundamental vector fields L a and Lo in Eq.(Al9) 
one can now introduce the fundamental horizontal vector 
fields H a and H on O(M) for the extended stochastic 
differential system Eq.(B2): 



Ha ~ ea dx* 



H = A i (r,X)d i -T i ml e' a (r)A' 



del, 



(B3) 
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and the operator Aq(m) f° r the stochastic process in the 
orthonormal frame bundle is given by 



A-o(m) = H a H a + Hq. (B4) 

a=l 

Aq(m) is the horizontal lift of the diffusion generator 
Am on the manifold to the orthonormal frame bundle. 
Obviously, the projection of a function in 0(M) to M 
with /(r) = /(x, 0), r = (x 1 , e*) satisfy the relation 

A OC M)/(r) = A M /(x) (B5) 

where A M = § E'=i.(4^4^) + A% = (f A + 
and Am = g^didj + g^T^dk is the Laplace-Beltrami op- 
erator. The generalized Kolmogorov backward equation 
on a Ricmannian manifold is obtained by 

|-u(r, x) = (^A M + A%)u(r, x) (B6) 

As shown in the appendix A the generalized Fokker- 
Planck equation is given by the adjoint of the diffusion 
generator A* (which includes the volume element ^/g, 
g = det{<7ij}). Since the Laplace-Beltrami operator is 
sclf-adjoint Am = A* M the generalized Fokker-Plack 
equation on a Ricmannian manifold takes the form: 

d$ D 

— = -div x LA$) + -A M $, (B7) 

where div x (A$>) = g~^ di(g^A l $) is the divergence op- 
erator in the Riemannian manifold, $ = $(x, t | y,0) 
is the transition probability with the initial condition 
$(x, | y,0) = <5(x-y) and adequate boundary condi- 
tions at infinity. The probability density </?(x, r) is de- 
termined by the same equation with the initial condition 
^(x,r = 0) = ^°(x). 
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